B-Spline Based Numerical Algorithm for Singularly Perturbed Problem of Fourth Order
نویسندگان
چکیده
In the present paper, a numerical algorithm using fifth degree quintic B-spline for fourth order singular perturbation problem has been developed. The most of the numerical methods used for higher order singularly perturbed boundary value problems transform the problems into equivalent system of first and/or second order differential equations. However, in the present method, fifth degree B-spline is applied directly to the problem without transforming the problem into an equivalent system. The method uses values of fifth degree B-spline function and its derivatives up to the order four at nodal points. Resulting system of equations is solved to get the required quintic B-spline solution. Since perturbed problems contain boundary layers, the strategy of fitted mesh is used which assigns more mesh points in the boundary layer regions. The algorithm is tested on two problems to demonstrate the practical usefulness and superiority of the approach.
منابع مشابه
Numerical method for singularly perturbed fourth order ordinary differential equations of convection-diffusion type
In this paper, we have proposed a numerical method for singularly perturbed fourth order ordinary differential equations of convection-diffusion type. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference method. In order to get a numerical solution for the derivative of the solution, the given interval is divided in...
متن کاملSolution of a Fourth Order Singularly Perturbed Boundary Value Problem Using Quintic Spline
Singularly perturbed boundary value problem can be solved using various techniques. The solution of the following fourth order self adjoint singularly perturbed boundary value problem is approximated using quintic spline Ly = − y(4) + p(x)y = f(x), p(x) ≥ p > 0, y(a) = α0, y(b) = α1, y(1)(a) = α2, y(1)(b) = α3, } or Ly = − y(4) + p(x)y = f(x), p(x) ≥ p > 0, y(a) = α0, y(b) = α1, y(2)(a) = α4, y...
متن کاملA Parameter Uniform Numerical Scheme for Singularly Perturbed Differential-difference Equations with Mixed Shifts
In this paper, we consider a second-order singularly perturbed differential-difference equations with mixed delay and advance parameters. At first, we approximate the model problem by an upwind finite difference scheme on a Shishkin mesh. We know that the upwind scheme is stable and its solution is oscillation free, but it gives lower order of accuracy. So, to increase the convergence, we propo...
متن کاملAn efficient numerical method for singularly perturbed second order ordinary differential equation
In this paper an exponentially fitted finite difference method is presented for solving singularly perturbed two-point boundary value problems with the boundary layer. A fitting factor is introduced and the model equation is discretized by a finite difference scheme on an uniform mesh. Thomas algorithm is used to solve the tri-diagonal system. The stability of the algorithm is investigated. It ...
متن کاملA parametric spline method for second-order singularly perturbed boundary-value problem
A numerical method based on parametric spline with adaptive parameter is given for the secondorder singularly perturbed two-point boundary value problems of the form 0 1 ( ) ( ) ( ); ( ) ; ( ) y p x y q x y r x y a y b The derived method is second-order and fourth-order convergence depending on the choice of the two parameters and . Error analysis of a method is briefly...
متن کامل